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^ ! Abstract 

Qh' We study the behavior at infinity, with respect to the space variable, of solutions 

■^ . to the magnetohydrodynamics equations in W^. We prove that if the initial magnetic 

fT^ I field decays sufficiently fast, then the plasma flow behaves as a solution of the free 

nonstationnary Navier-Stokes equations when |x| — > +00, and that the magnetic field 
will govern the decay of the plasma, if it is poorly localized at the beginning of the evo- 
Ph i lution. Our main tools are boundedness criteria for convolution operators in weighted 

spaces. 
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^ '■ 1 Introduction 

cn I The magnetohydrodynamics equations are a well-known model in plasma physics, describing 

Q ■ the interactions between a magnetic field and a fiuid made of moving electrically charged 

\0 • particles. A common example of an application of this model is the design of tokamaks: 

^ ! the purpose of these machines is to confine a plasma in a region, with a density and a 

^ I temperature large enough to entertain thermonuclear fusion reactions. This can be achieved, 



at least during a small time interval, by applying strong magnetic fields. We refer to [12] 
for other applications of this model, in particular to the study of the dynamics of the solar 

K^ ■ corona. 

K^ , In non-dimensional form, the magnetohydrodynamics equations can be written in the 

H [ following way: 

^ ^ + {u-V)u- S{B .V)B + v(p + ^\B\A = ^Au 

f)R 1 

-- + («■ V)B - {B ■ V)u = -—AB (MHD) 

or -ti-m 

div u = div B = 
^ u{0) = Mo and 5(0) = Bq. 

Here the unknowns are the velocity field u of the fiuid, the pressure p and the magnetic 
field B, all defined in M'^ {d > 2). The positive constants Re and Rm are respectively the 
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Reynolds number and the magnetic Reynolds number; moreover S = M'^/{ReRm), where M 
is the Hartman number. After rescaling u and B, we can assume that S = Re = I. With 
minor loss of generality, from now on we shall also assume that Rm = 1. All the results 
however remain valid in the general case with simple modifications in the constants. 

In the particular case 5 = 0, the system (MHD) reduces to the celebrated Navier-Stokes 
equations. Just as in this particular case, global weak solutions to (MHD) do exist, but 
their unicity, as well their smoothness in the case of smooth data, remains an open problem 
for d > 3. Partial regularity results, which provide bounds of the Hausdorff dimension 
of the possible singular set of weak solutions, have been obtained in [7]. Constantin and 
Fefferman's theory [5] relating the regularity of the fiow to the directions of the vorticity 
has been extended to magnetohydrodynamics in [8]. A construction of forward selfsimilar 
solutions is given in [9], where the nonexistence of backward selfsimilar solutions is also 
discussed. Moreover, the asymptotic behavior of the solutions for t -^ +cxd is quite well 
understood: for example, [13] provides the optimal decay rates of the L^ norm of u and B 
for a large class of fiows. 

On the other hand, nothing seems to have been done to study the decay of solutions 
of (MHD) with respect to the space variable. In this paper, motivated by recent results 
obtained by several authors for the Navier-Stokes equations (see, e.g., [1], [2], [6], [11] and 
[14]), we would like to describe in which way the presence of the magnetic field affects the 
spatial localization of the velocity field. 

Definitions and notations. We start by introducing the notion of decay rate at infinity 
in a weak sense, which generalizes the usual notion of pointwise decay rate in the framework 
of locally square integrable functions. A simple motivation is that the L^^^ regularity is the 
minimal one for which the system (MHD) makes sense. 

1. Let / G L^Qp(M'^). We define the L^ decay rate as \x\ —^ +oo of /, as 

ri{f) = snp\rieR; lira R^'^f \f{Rx)\^dx = o\. (1.1) 

I R^+oo Jl<\x\<2 ) 

li rj = ri{f) is finite then we will write / ~ |a;|~'' when |a;| -^ +oo. On the other hand, 
when we write / = 0(|a:|~'') when |x| -^ +oo, we mean that ?](/) > rj. Of course, 
any measurable function such that |/(a;)| < C(l + |a;|)~'' satisfies / = (9(|a;|~'') when 

\x\ -^ +CXD. 

2. For a e [1, +oo] and a G M, the space L'^iM.'^) is the Banach space normed by 

ll/lUg = U l/(^)r(l + 1^1)"" dx] if 1 < o < +00 (1.2a) 

and, if a = +oo, by 

LOO = ess sup |/(x)|(l + |x|)". (1.2b) 



From the localization point of view the two spaces LJ^(R'^) and L^(M'^) must be consid- 
ered as equivalent, when 

d ^ d 

a + - = /?+-. 
a 
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Indeed, if / G L''^{R'^) and a > 2, then / = C> ( | a; | -("+'='/")) when |x| -^ +00. Holder 
inequality implies that 



L^a C L^ (1.3) 

whenever a + d/a > P + d/b and a > b. It also implies that 

r7(/)=sup|a + ^; a>2 and /gl4 (1.4) 

for any / G Ll^iM.''). 
We shall use the following additional notations : 

3. If A and B are two expressions containing a parameter a, then when we write 

A<B-ea, 

we mean that A < B ifa = and A < B iia j^ 0. We shall also often write expressions 
of the form A < B — ei/a meaning that the inequality must be strict for finite a and 
can be large when a = +00. 

4. The positive part of a real number will be denoted by (■)"'' = max{-, 0}. 

Main results. We are concerned with the persistence problem of the spatial localization 
of the magnetic and the velocity fields. Our main results (Theorem 1.1 and 1.3 below) aim 
to answer the following questions. Consider a localization condition like 

{uo,Bo)eL2(m')xLl{R'). (1.5) 

Will the unique solution of (MHD) preserve such a condition in some future time interval ? 
Depending on the parameters, the answer can be positive or negative. In case of a negative 
answer, can we still ensure that the spatial localization of the solution is conserved in the 
weak sense ? In other words, we would like to know whether 



-00. 



u{t) = 0{\x\'^^''+'^/P"^) and B{t) = 0{\x\-^^'+'^/p'^) when |x| - 

Again, this condition may be conserved, or instantaneously break down. 
We will prove the following: 

Theorem 1.1 Let uq G I^^°(M'^), Bq G Ll\(M'^) be two divergence-free vector fields in R'^ 
(d >2). Assume that 

^0 > f ^1 > ^ ^ 

and < (1.6a) 

d < po < +00 \ d < pi < +00. 

Let us also assume that 

f^ + e^^ < r^o < min{rf+ 1; 2?7i -(5}, (1.6b) 

with r^o = f^o + d/po, % = "^1 + d/pi and 5 = i— — l\ . Finally, define p^ = min{po ; f — ^5} . 



Then there exists T > and a unique mild solution [u, B) of (MHD) in C( [0, T] ; L^o xLP'). 
This solution satisfies 

«(t) = C(|xr*) and B{t) = 0{\x\-"') when \x\ -^ +00. (1.7) 

ff d = 2, the time T can be arbitrarily large. 

Moreover, if {uq, Bq) also belongs to L~° x L~^ , with the corresponding indices satisfying 



sup 

te[o,T] 


\u{t. 


x) 


L_2 


0{x 


-(d+l+e)\ 


sup 

te[o,T] 


Bit,. 


x) 


L_2 


0{x 


-id+l+e)/2\ 



assumptions (1.6), then the lifetimes in L^o x L^^ and L^" x L^^ agree and both maximal 
solutions are actually the same one. 

Next we discuss the optimality of the above restrictions. Such restrictions are of two 
kinds: there are a few conditions related to the well-posedness of the system, and a condition 
(namely, the upper bound for rjo in (1.6b)) which is related to the spatial localization of the 
solution. Here, we will only focus on this condition. The following theorem implies that the 
restriction tjq < d+1 is sharp. We expect that the other restriction is also sharp, or at least 
that ?7o < 2?7i for stable weak solutions. But we were not able to prove such a result. 

Theorem 1.2 Let {u,B) e C([0,T]; ^^(M^) x L'^{R'^)) a solution to (MHD) such that 

(1.8a) 

and snp \B{t,x)\ = 0{\x\-'^''+'+'>^'') (1.8b) 

te[o,T] 

for some e > 0. Then, for all t G [0,T], there exists a constant C{t) > such that the 
components ofu{t) and B(t) satisfy the following integral identity : 

f {u^u'' - B^B^){t, x) dx = 5j,k C{t), {j,k = l,..., d) (1.9) 

with 6j^k = 1 ifj = k and 6j^k = otherwise. 

By Theorem 1.3 below, condition (1.8b) will be fulfilled as soon as mq and Bq belong to 
Ll{W'-), with p> d and f? + - = (rf + 1 + e)/2, for some e > 0. This means that if we start 
with a well localized initial datum {uq, Bq), but such that (1.9) does not hold for t = 0, then 
condition (1.8a) must brake down. 

On the other hand, the integral identities (1.9) are obviously unstable. Neverthless, in 
section 5 we shall see that a class of exceptional solutions satisfying (1.9) does exist. Inside 
this class, one can exhibit solutions such that u decays much faster than in the generic case. 

Physical interpretation of Theorem 1.1. This theorem reinforces mathematically some 
facts that can be observed in the applications. Three conclusions can be drawn: 

1. Any spatial localization assumption on the magnetic field will be conserved by the 
flow. Indeed, the L^ decay rate 771 can be arbitrarily large. The spatial localization of 
the velocity field is also conserved, but there are some limitations to this property. 



d+1 




\vi<id + l)/2 
[pi > 2d 
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5/d 



1/d 2/pi 



'm < {d+l + 5)/2 
^d <pi <2d 



2rii 
d+l 




1/d p 

jm > {d + i + 6)/2 

\pi >d 



Fig.l The figures show the admissible 
values for (po?^o) allowing (1.7) to hold, 
once (pi, -i^i) is given {all gray regions). 

Above : Slowly decaying magnetic field. 
The results depends slightly on the regu- 
larity of B through 6 = ( ^ - 1 

Down-Left : Fast decaying magnetic 
field. The velocity field behaves at infin- 
ity as the solution of Navier-Stokes equa- 
tions with the same initial datum Uq (see 
[14]). 

The dark gray regions correspond to ini- 
tial data for wich we will prove in addi- 
tion that u e L°°{[0,T];L2). The dash- 
dotted lines illustrate the barriers used in 
the proof of §4.3. 



2. For poorly localized magnetic fields (namely rji < {d + 1 + S)/2), the behavior of u 
when \x\ -^ +oo is governed by the decay of the magnetic field. As < 5 < 1 in (1.6b), 
the maximal L^ decay rate of u that can be conserved by the fiow exceeds 2?7i — 1. 
When pi > 2d, one has 6 = and this rate is improved up to twice that of Bq. The 
pathological lower bound on rjo disappears too. Roughly speaking, requiring pi to be 
larger (for a given L^ decay rate rji = {}i + d/pi of the magnetic field) means that the 



behavior at infinity of Bq is closer and closer to that of a function that decays as |x| ''^ 
in the usual pointwise sense. 

3. For sufficiently fast decaying magnetic fields, the decay of u is not affected by B, but 
is provided by the fundamental laws of hydrodynamics. The reason is the following: 
for magnetic fields such that rji > {d + 1 + S)/2, our limitations on the L^ decay rate 
at infinity of the velocity field (1.6b) boil down to the only restriction rjo ^ d+1. This 
is exactly the same restriction that appears for the Navier-Stokes equations. Indeed, 
we know from F. Vigneron's result [14] that the mild solution of the Navier-Stokes 
equations remains in L^° (M'') if the initial velocity belongs to such space and 

A + d/pQ <d+l -£i/po- 

This condition in known to be sharp. One may notice however that, thanks to (1.4), 
the equality case is possible even if po < +oo, provided that stability is asserted as in 

(1.7). 

A more physical explanation for the above conclusions is the following^. The induction 
equation means that the magnetic field lines are transported by the fiow while simultaneously 
undergoing resistive diffusion. This transport-diffusion process guarantees that, where the 
velocity vanishes, the magnetic field will not spatially spread out during small time intervals, 
since the mechanism of diffusion is quite slow. As for the fiuid fiow, the magnetic field acts 
upon it only through the Lorentz force: whenever this disappears the velocity acts in a purely 
Navier-Stokes way; thus, the spatial spreading of the initial velocity is essentially governed 
by the competition between diffusion, whose effect is important only for large time, and 
incompressibility, that immediately prevents the fiow from remaining too localized. 

Stability in weighted spaces. Conclusion (1.7) does not mean that 

(«,5)GL-([0,T];L«;xL^,;). 

Actually, we do not know if this property holds when uq G LJJ° and (po, ^o) is in the light-gray 
regions of Fig.l. However, if (pc^o) is in a dark-gray region, then such property does hold. 
This is essentially the statement of our next theorem. It extends to the case of non- vanishing 
magnetic fields, the result established in [14] for the Navier-Stokes equations. 

Theorem 1.3 Let uq e L^iM'^), Bq e I/5'^(M'^) be two divergence-free vector Relds in W^ 
(d > 2). Assume that i^o, ^i > 0, d < po < +oo and 

2 11 

— <— + -. (1.10a) 

Pi Po d 

Then there exist T > (if d = 2, one may take T = +ooJ and a unique mild solution of 
(MHD) 

(m,5) gC([0,T];L«' xLPi). (1.10b) 

If, in addition, the decay rates of Uq and Bq defined by rjo = {^q + d/po and r/i = ^i + d/pi 
satisfy 

f , d 2d] , , 

?7o < mm<^ d+1- ei/p„ ; 2r]i - £2^1-^0 ; Sr^i H }, (1.11a) 

I Po Pi J 

^This explanation was suggested to us by the Referee. 



then we have more precisely 

iu,B)eC{[0,T];L2xLl\). (1.11b) 

Moreover, if {uq, Bq) also belongs to LF~ x L~\ with new indices again satisfying (1.10a) 

and (1.11a), then the lifetimes in L^" x L^^ and L~° x L~^ are the same and both maximal 
solutions agree. 

The assumption (1.10a) is not really related to spatial localization problems, but rather 
to well-posedness issues of the equations, and in particular, to the invariance of the equation 
under the natural scaling 

uxit, x) = \u{\H, \x), Exit, x) = \B{\H, \x) (A > 0). 

We expect that Theorem 1.3 remains true in limit cases p=d, or — = — + 2 (with several 
modifications in the proof). We did not treat these limit cases since they would require Kato's 
two-norm approach for proving the boundedness of the operators involved, as described in 
[3, chap. 3] or [4] for the Navier-Stokes equations. The proof would be more complicated, 
without providing any substantial clarification of the spatial localization problem. 

Let us also observe that one could replace the weights (1 + |a;|)'' with homogeneous 
weights. But in this case the conditions to be imposed on the parameters would be much 
more restrictive, e.g. 

^ + -<1. 
P 

Again, this would not help to understand the spatial localization of the fields. 

Main methods and organization of the paper. We shall first prove Theorem 1.3 and 
later deduce Theorem 1.1 as a corollary of the natural embedding (1.3) between weighted 
spaces. The idea consists in observing that the assumptions (1.6), together with the inclusion 
(1.3), ensure that the initial datum belongs to the product of two larger Lebesgue spaces, in 
which we can prove the existence and uniqueness of a mild solution. 

Our proof of Theorem 1.3 consists in applying the contraction mapping principle to the 
integral form of (MHD), in a suitable ball of the space C([0,T],L^° x L^\). This is why we 
refer to {u, B) as a mild solution. The only difficulty is establishing the bicontinuity of the 
bilinear operator involved. 

For small values of ?7o, the bicontinuity would be a straightforward consequence of the 
well-known Young convolution inequality in weighted Lebesgue spaces (recalled in [14, §2.2]). 
But this argument does not go through when r/o is close to the upper bound of (1.11a), since 
the kernel of the operator governing the evolution of the velocity field decays too slowly at 
infinity. In this case, the proof requires more careful estimates. The main one is given by 
Proposition 3.1 below. 

Several generalizations of the weighted convolution inequalities are known (see, e.g., the 
recent boundedness criterion for asymmetric kernel operators [14, §2.3], which applies to 
Navier-Stokes) . However, we could not deduce the bicontinuity of the bilinear operator by 
applying directly any known inequality, unless we put additional artificial restrictions on the 
parameters. 

The main issue with the spatial localization of magnetohydrodynamics fields is that the 
system cannot be treated as a scalar equation. When dealing with the Navier-Stokes system. 



one may often reduce the problem to a single equation, because all the components of the 
kernels of the Navier-Stokes operators satisfy the same estimates. This is no longer true 
for (MHD). In the following, we shall derive sharp bounds for the magnetohydrodynamics 
kernels and take advantage of the fact that a few components decay much faster than the 
others. 

This paper is organized as follows. Section 2 contains some generalities on magnetohy- 
drodynamics. In Section 3 we study the boundedness of convolution operators in weighted 
spaces. We use these results in Section 4, proving first the local existence of a unique solution 
in weighted spaces (1.11b), then the fact that lifetimes do not depend on the choice of the 
indices. Then we deduce Theorem 1.1 as a corollary. 

Theorem 1.2 will be proved in Section 5, using a Fourier transform method developed in 
[2]. Section 5 also contains the description of a method for obtaining special solutions, such 
that the velocity field is more localized than in (1.6b). Those solutions are however unstable. 

Remark 1.4 When we deal with the space C([0,T];L^° x L^^^), with po = +oo or pi = 
+00, the continuity at t = must be understood in the weak sense, as is usually done in 
nonseparable spaces. 

2 The integral form of the equations 

Let P be the Leray-Hopf projector onto the divergence-free vector field, defined by 

P/ = /-VA'i(div/). 

Applying P to the first equation of (MHD) and then the Duhamel formula, we obtain the 
integral equations 

u{t) = e'^Uo- / e^'-'^^Fdiv{u^u-B®B){s)ds 
Jo 

B{t) = e'^Bo - [ e(*-^)^ dw{u (» B - B ^ u){s) ds ^^^^ 

Jo 

divMo = div_Bo = 

where e*"^ is the heat semigroup (recall that the Reynolds numbers and the Hartman numbers 
have been set equal to 1). The semigroup method that we use in this paper to solve (IE) 
provides mild solutions of (MHD) that are in fact smooth for strictly positive t. 

We denote respectively by F^hit, x) and C^f^it, x) (j, h,k = 1, . . . ,d) the components of 
the kernels of the matricial operators e*^PV and e'^^V. Thus, 

^>(e,t) = e-*i«i'a(5,>-e,e.ier'). (2.1) 

This expression of the symbol allows us to see that 

F{t,x) = t-^'^'y'^x/Vt), 

with |$(x)| <C(l + |x|)-(^+^). 



< 



This low decay rate of $ is due to the fact that F(t, ■) ^ Lj(M'^) ; otherwise F(t, ■) would be 
a C^ function on W^. On the other hand, 

, V / ;, 2.2b 

with ^ e S{R'^) (the Schwartz class). 
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Let US introduce the bilinear operators on M*^- vector fields U and B whose k^^ component 

l]'{f,g){t,x) = Y^f Fl^{t-s)*{r®g^){s)ds 
j,h ■'^ 

M'{f,g){t,x) = Y.fGl^{t-s)*{r®g^){s)ds, 
and the bilinear operator V = (Vi, V2) on M^'^-vector fields v = (^1,^2) defined by 

¥i{v,w) = V{vuWi)-U{v2,W2) 
Y2iv,w) = M{VUW2)-M{V2,W,). 

Here and below, for v G M?'^, we denote by Vi the first d components and by V2 the last d 
components. 

With these notations and setting v = {u,B), vq = {uo,Bq), the system (IE) can be 
rewritten as 

v = e*^vo-Y{v,v). (2.3) 

As it is well known (we refer, e.g., to [3, Lemma 1.2.6]), if X is a Banach space, then for 
solving an equation like (2.3) one just needs to check that 

e'\o^C{[0,T];X) (2.4a) 

and 

V : C([0,T];X)xC([0,T];X)^C([0,T];X), (2.4b) 

with the operator norm of V tending to as T ^ 0. Then the existence of a solution 
V G C{[0,T];X) is ensured, at least for T > small enough. 

In order to prove Theorem 1.3 we shall take X = L^° x L^^^. In this setting, condition 
(2.4a), the unicity and the continuity of the solution with respect to the time variable are all 
straightforward. Therefore, our attention will now be exclusively devoted to the more subtle 
problem of the bicontinuity of V in L°°([0, T]; L^° x Lf^). 

We need three estimates, namely 

||U(M,M)(t)||iPoA. < '^t||m|Ic([o,t],lpo.^o) (2-5a) 

||U(5,5)(t)|L™,.„ < Cr||i?||^([o,T],L.i.^x) (2.5b) 

[u,B){t)\\i^Pi,«-^ < C'r||M||c([o,T],LPo.''o)||-B||c([o,r],LPi'''i) (2.5c) 



for all < t < T and some constant Ct such that Ct — > as T ^ 0. These bounds will not 
rely on the specific structure of the operators U and B, but only on the decay properties of 
their respective kernels: 



|F(t,x)|<C(v^+|x|)-('^+^) 



■N~d-l 



\G{t,x)\<CNVt (v^+|x|) 

for all iV > 0. 

We start by observing that by Holder inequality, 



-N 



(2.6) 



\u®u\ ,po/2 


< \u 


|2 

■Bq 


B®B\ ^pi/2 


< \B 


1 2 

^^1 


B Lh(po,pi) 


< \u 
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L'.' 



where -^^^ — y = — + — denotes the Holder exponent (the assumptions of Theorem 1.3 imply 
that po,pi > 2). Set A = \/t and 

r^(x) = (A + N)-^. (2.7) 

Then the only thing that we have to do to obtain (2.5a)-(2.5c) is to establish that for all 
0< A< 1: 

||rf^*/|L.o<CA- 11/11^.0/2, (2.8a) 



^2i9 







\\Tf'*fhf<CX^oiif\i (2.8b) 

and ||r^ * /IL.1 < CA'^^ll/ll^ec.o,,,) (2.8c) 

with an arbitrarily large A^ > and exponent (Tq, ctq, (Ti such that 

(Jo > -2, (Jo > -2, ai> -N + d- 1. (2.9) 

The constant C > has to be independent of A. Assumption (2.9) ensures that the integrals 

/ \\F{t-s)*{u0u){s)\\Ln,ds, I \\F{t-s)*{B®B){s)\yods 

Jo ''o Jo "0 

and 



\\G{t - s) * {u ® B) {s)hr., ds 

'0 ^1 

converge. 

3 Convolution estimates in weighted spaces 

The fundamental estimates (2.8a)-(2.8c) will be a simple consequence of the following propo- 
sition. 
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Proposition 3.1 Let a,p G [1; +00] and a,^ > 0. For any real numbers A > and N > 1 
let us set 

T^ix) = iX + \x\r^. 

Let also f e L''^(R'^) and N > d. 

1. Then T^ * f e /.^(M'^), provided that 

■& < a and ■& -\ — < min< N — ei/p ; a H Sa--& (■ (3.1) 

Moreover, if N j^ d{l + - — ^) , then there exists C > such that 

||rf*/||,.<CA-^(l + A)^||/|U,. (3.2) 

2. If one assumes in addition that 

111 , , 

-<- + -, 3.3 

a p a 

then there exists e > and two constants C,m > such that 

lirf * /Iks < cx-'^^'-'^^ii + xrWfllK- (3.4) 

When N = rf(l + i-i), the bounds (3.2) and (3.4) hold with an additional factor (1 + | logA|) 
in the right-hand sides. In (3.2) and (3.4) the constant C may depend on d, a, a, N and d, 
but it does not depend on A or f. 

Remark 3.2 We shall see in the proof that we can take 

(d d N-d + 1 

e = min < hi; 

{pa 2 

m = max \n -d + l-2e; -N + d I hi 

Proof. We start by observing that by Holder's inequality, 

<C||/|U. ii-<-<mm\l;- + ^-ea 

a q [ad 



Next we have 

{i + \x\r\T^*f{x)\< 

with the following definitions : 



(3.5) 



(1 + \x\f = I^^xix) + J.ff,xix) + K^,x{x), 



T^,{x-y)\f{y)\dy 



hAx)=( / T';!{x-y)\f{y)\dy Ul + \x\r, 

V"'|y|>|2:|/2 / 

Axix) = ( f T^ix-y)\fiy)\dy)il + \x\fln^o,i)ix), 

V-^ls/|<l^l/2 / 



K^Ax)=[ / T';!{x-y)\f{y)\dy\{l + \x\riBioM^)- 

\J\y\<\x\/2 J 

Here and below, 5(0, 1) denotes the unit ball and \e is the indicator function of a set 

E C M^. 
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The bound for K.g^\. Since \y\ < \x\/2, we have 

(A+|x-i/|)-^<2^(A+|x|)-^. 
Hence, using (3.5) with i^ = 1 - i = (i _ | - i + 5^) + , 

J\y\<^4 
<C(A + N)-(^-^)||/|U, ||l^(o,H/2)IL,' 

< C(A+ |x|)-(^-'') k|f'-("+«)+^"l^||/||Lg. 



As \x\ > 1, it follows that H-ft'ij.AllLP < C* 11/11 l°, uniformly for A > 0, provided that 



d — { a + 



P 



Since N > d, this condition is weaker than (3.1). 



£i/p- 



The bound for J^^x. Using (3.5) again, but with q = a, gives us 



0<X&A^)<ClB(o,i)ix)i\+\x\r'' / \fiy)\dy 

< Cl^(o,i)(:r) (A+ |x|)-^H'^(i-i/'^) ll/IU., 



(3.6) 



whence 












1 J'd,\ \lv < C 


A-^^ / \X ^(l-l/'^) dx + 1{A<1} 


4 


|^|-iVp+dp(l-l/a)^^ 

x|<l 


i/p 


I/Il- 


Thus, for all i? > and p G [1, +00], we have 








J,^\l.<c[i + \-^^'^'^--) i/u. 




/ 1 1\ 

ifA^^rf 1 + , 

\ p aJ 


(3.7a) 


and J^ 


aIl. <C(1+ logAI) l/U. 




( 1 l^ 

\i N = d 1 + 

V p ay 




(3.7b) 



Note that || J^^aIIlp is bounded by the right-hand side of (3.2). Moreover, if ^ < - + ^, then 
II'^'iJ.aIIlp is also bounded by the right-hand side of (3.4), provided that < e < (i(- — ^ + ^)- 

The bound for U^x- Set F[x) = (1 + \x\Y |/(x)|, so that F G /."(R'^) and 

0</.,A(:r)<C(l + |x|)-("-^) /" Yl{x-y)F{y)dy. 



But rf e L^(M'^) for all 6 G [1; +00] and /3 > such that /? + f < A^ - ^i/fe. Moreover, one 
has 

l|r^llL^<CA-"^+f(l + A)^. (3.8) 

The remaining part of the proof of Proposition 3.1 relies on the following lemma. 
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Lemma 3.3 Let a,b,p e [l;+oo] and a,p,^> 0. For f G L''^{R'^), g E L^lM'^), define 

I^x) = {I + \x\)-^''-'^^ F * g{x) 
with F{x) = (1 + |a;|)" |/(ic)|. If there exists s G [1, +oo] such that: 
( ^<a 



d \ d f d 

— < mill < — ; \ a -\ — 
s \ a \ a 



. d\ ,/ 1 

V + -\ - Ea^-O ; rf I 1 - - 



(3.9a) 



d \ d d 
- > max < ; 

s \ 0. p 



d- [P+-] +ep 



then U G L^iW^) and 



(3.9b) 



Proof. According to (3.5), we have g E L^ {W^) for all s' G [1; +oo] such that 

11 r 1 /? 

Let J + ^ = 1. We now use that ^ — -^ > 0. The Young exponent Y(a, s') of a and s' is well 

_ _ 1 _ 1 

') a s- 



defined by yiIs') ^ \~\- Moreover, one has F * g e L^'-^'^'\W'-), i.e. 






Since ^ < a, (3.5) implies that /^ G U'iW'-) for all p such that 



111 r 1 1 a-'d 

< - < mm <l] \ e„_^ j> , 

asp y a s d 



and (3.9b) is satisfied. 



D 



Let us now come back to the proof of Proposition 3.1. We are going to apply the lemma 
with g = T^ , I^ = I^ Xi b = +oo and P = N. 



If - < -, then we further choose s = +oo and conditions (3.9a) boil down (recall that 
A^ > d) to the only restriction ■& + -<a + - — ea--d- 



- If - > -, then we choose - = - — -. In this case conditions (3.9a) boil down to ^ < a. 

ap' sap V/ — 

The first part of Proposition 3.1 now follows from the bounds obtained for /^ a, J'd,\ and 

To prove (3.4), we fix e such that < e < ^-^+'^ . Then we apply Lemma 3.3 again with 
g = T^ and /^ = /^ a, but with b and P defined by 



d 



d-l + e, and p = N-d+l-2e. 



By (3.8), one has Pf G L''JR'^) with ||Pf l^o < A" 



^0(A)and0GL-([O;+oo)). 
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As before, 

^ if ^ < -, then we choose s = +00 in (3.9a) and Lemma 3.3 implies that 

-N+d~l+e 



/.,a|L^<A-^^+^-^+^0(A)1|/1U., (3.10) 



"/3 



provided that 'd+-<a + ^~ €a-^. 
- If i > i, then 1 = 1-1 leads again to (3.10), provided that ^ < a and 1 < 1 + 1 - -1. 

a p^ sap o \ /I L- — (J — pad 

The proof of Proposition 3.1 is now complete. D 

4 End of the proof of Theorems 1.1 and 1.3 

4.1 Existence of a unique mild solution in weighted spaces 

We are now in a position to prove Theorem 1.3. 

Under the assumptions of Theorem 1.3, one applies (3.4) with N = d + 1 and with 
e = l — ^ore = l— (— — -) respectively; assumption (3.3) is ensured by (1.10a). This 
proves (2.8a) and (2.8b) with 

d . , f2d d^^ 



(To = — 1 and 

Pq " \Pi Po 

A new application of (3.4) with any A^ such that A^ > raa,x{d + 1 ; ^i + — } + Si/p-^ and 
e = 1 - -^ yields (2.8c) with ai = -N + d - d/po. 

With the preceeding values of o-QjCTq and ai, the assumption (1.10a) implies (2.9). As 
indicated in section 2, this yields (2.4b) and ensures that the operator norm of V tends to 
zero as a power of T, when T ^ : 

< C max i Ti+^ ; T'^"^ ; T'+hi-^+^-d-D 



C{[0,T];X) 

This ensures finally the conclusions (1.10b) and (1.11b) of Theorem 1.3. 

More precisely, our argument proves that under the assumptions of Theorem 1.3, the 
maximal lifetime T* of the mild solution in X = L^' x L^^^ satisfies 

{-2/(1--^) -2/('i-r^-^i^'i 1 

l;||(«o,5o)|lx ''^ ;||(«o,5o)|lx^ ^" ™^ ^[, (4.1) 

with a constant c > 0, depending on all the parameters, but not on Uq or on Bq. 

4.2 Comparison of lifetimes in Theorem 1.3 

It only remains to establish that lifetimes are independent of the admissible pairs of indices 
chosen to construct the solution. 
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Proposition 4.1 Let uq G L^i^"^) n L|«(M'^) and Bo E Ll\{R'^). Set r/o = ^o + d/po, 
ifo = ^0 + d/po and rji = {}i + d/pi. Assume that d> 2 and 

d <Po,Po < +00 

2 r 1 1 1 1 

— < mill < h - ; Tc- + - 

Pi Ipo " Po « 

r , d 2d] (4.2) 

?7o < mm <^ d + 1 - el/p^, ; 2?7i - £2^1-^1 ; Sr^i H 

I Po Pi 

^ f , rf 2rf 

r^o < mm <^ d + 1 - ei/p- ; 2r7i - e^„ _g ; 2r7i + — 

[ ^ " Po Pi 

Let T* and T be the Metina.es of the solution {u, B) of (mhd) emanating from (mq, -Bo) in the 
respective weighted spaces, i.e. 

T*=sup{T>0 s.t. {u,B)eC{[0,T];L2xLl\)}, 
f = sup I T > s.t. {u, B) e C([0, T]; L| x Ll\) } . 

Then f = T*. 

Proof. The structure of the proof is similar to that of [14]. Let us assume that we have, for 
example, T < T*. Unicity of mild solutions ensures that they agree on [0,T[. We are going 
to prove that 



sup ||«(t)||^Po + l|5(t)||^Pi < +00. 
iG[0,T[ V ''0 V 

Then (4.1) would imply that the mild solution (u, B) in L~° x L^^^ could be extended beyond 

T, and that would contradict the definition of T. 

First of all, let us recall (see, e.g., [14, §2.2]) that there exists a constant Cq > depending 
only on d and ■(?, such that 

sup VvW^v. < Co (1 + fy^" \\v\y. . (4.3) 

re[0,T] 

In the following, we set A = Co (1 + fy^/'^. 

Note also that we can obviously assume that m ^ in [0, T]. 

The bound for B. By the second of the integral equations (IE), one has for < s < t < T : 

B{t) = e(*-")^5(s) - I G{t - t) * {u® B - B ® u) {t) dr. 

Proposition 3.1 applied to the upper bound of G given by (2.6), with e = 1 — ^ in (3.4), 
yields 

yT<t<f, ||G(t-r)*(M®5)(r)||^Pi <ir(t-r)-'^||(M®5)(r)|Le(p„.p,) 
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where a = ^{1 -\ ) and i^ is a constant, possibly depending on T* and all the parameters 

contained in (4.2), but not on T. Note that a < 1. Thus, for all t G [0;T], 

||5(t)|L.. <A\\Bis)h.r +K ^^~'^ sup \Ht)\\^.o ■ sup ||5(r)|Lpi . (4.4) 

"i "i 1-0- re[s,t] "» TG[s,t] ''i 

Now let (Tn)n>o be the increasing sequence defined by 

Tn = nA with A= ( sup ||M(r)||^po 

re[o,f] 

and iV G N such that T^ <f < T^+i. FovO <n<N, let /„ be the interval [T„, T„+i] n [0, f[ 
and 

M„ = sup||5(r)||^Pi. 

Applying (4.4) with s = Tn and t G /„ for n = 0, . . . , A^, we get 

Mo<2A\\Bo\\rPi and M„ < 2AM„_i (1 < n < iV), 
whence 

SUD ll^miL 

"''l 0<n<Af 



sup ||5(t)||iPi = max M„ < (2A)^+^||5o||iPi 



te[o,T[ 
Finally, this leads to : 

sup \\B{t)\\^.. <C||5o|L.iexp( fl+f sup ^(s)!!^!^^^) (l + ^i log(l + f ) 

(4.5) 
The right-hand side is finite because we assumed T < T*. 

The bound for u. For < s <t < T, one has 

u{t) = e^'~'^^u{s) - [ F{t-T)*{u0 u){t) dT+ [ F{t - t) * {B ® B){t) dr. 



Proposition 3.1, applied this time to the upper bound of F given by (2.6), yields 
||w(t)||^po <A\\u{s)\\^f,a +A'— sup ||M(r)||^po ■ sup ||M(r)||^po 

+ K ^^~ "V I sup ||5(r)||,.. 

with a = i(l + — ) and a = ^(1 + (— — — )+). Note that a is the same as before and 

that a < 1; K depends on T* and all the parameters, except T. The last term is uniformly 
bounded by 

KT^-^ ( 
L = . _~ sup^ \\B{t)\\lpi 

^ " \re[0,T[ ' 
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which is a finite constant because (4.5) holds. Define (T„)„>o and /„ as before. Let also 

Mn = sup ||M(r)||^po- 

Recall that N is the integer part of T/A. Then, for 1 < i < A^, one has 

Mo < 2A IImoII^po + 2L and M„ < 2A M„_i + 2L, 

hence 

sup ||M(t)|Lpi, = max M„ < (2A)^+^ II^oLpo + 2L [l + . . . + (2^)^-^ + (2A)^1 < +oo. 
te[o,f[ 5o o<"<^ % 

Combined with (4.1) and (4.5), this estimate ensures that T > T*. Exchanging the roles 
of T and T*, one finally obtains that T = T* . D 

An analogous result holds if we assume instead uq G L^° (M^) and Bq E L^\{W'-)r\L^^ C^"^), 
with obvious modifications in (4.2) : 

d < po < +00 



2 2 1 11 

max i — ;— >< \- -. 

Pi Pi) Po d 



(4.2') 



r]o < min {rf + 1 - ei/p^, ; 2r7i - £2^1-^0 ! 2r7i - €2-&,-^o} ■ 

f d 2d ^ d 2d] 

r]o < mm < 2r7i H ; 2r7i H -::- ^ • 

I Po Pi Po Pi J 

Theorem 1.3 is now established. 

4.3 The proof of Theorem 1.1 

Let po,Pi and ^oj^i such that (1.6a) and (1.6b) hold. 

If ^0 < 2^1, pq < d/6 — es and rjQ < d + 1 ~ £i/po, then (1.10a) and (1.11a) hold, and 
there is nothing more to prove since Theorem 1.3 already gives a stronger conclusion. 

In all the other cases and for any e > 0, our assumptions yield an embedding L^° C L^ 
such that Theorem 1.3 may be applied to 

and with 

d 

/^ + - = ^0 - e. 
Q 

It follows that u = 0(|x|-(*-^)) and B = 0{\x\-'^^) when \x\ -^ +00. Letting e -^ 0, this 
will conclude the proof of Theorem 1.1. 

Let us be more precise about the embedding L^° C L^. Actually, various choices are 
possible for (g, /i). We have chosen the indices that are represented on the interpolation 
diagram (see Fig. 1 p. 5) by a dash-dotted line. 
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If the magnetic field decays sufficiently fast, namely if ?7i > (rf + 1 + (5)/2, the only case 
not included in Theorem 1.3 is that of ?7o = c? + 1 with po finite. In this case, one may take 



(g,/^) = {Po,^o 



e . 



Let us now assume that rji < {d+ 1 + 6)/2 and, for the moment, that pi > 2d. Then the 
cases to be dealt with correspond either to ^o > 2^i or to rjQ = 2rii, or to both. 
-If^o >2^i, then 

- = ^0 - 2^1 + e and /x = 2^i 

Q Po 

are suitable, even if tjq = 2rji. 

- If -i^o < 2'/9i and r^o = 2rji, one may again choose (g, /i) = (po, ''^o — e)- 

Finally, if d < pi <2d and rji < (rf + 1 + 5)/2, one may use the following barrier : 

- = 1 — (1 — 5)k jj, = 2'i9i(l — k) and k, = 1 r— • 

q 2{r]i - 6) 

The proof of Theorem 1.1 is now complete. D 



5 Instantaneous spreading of rapidly decreasing fields 

This section is included for completeness and contains the proof of theorem 1.2, and some 
remarks about exceptional solutions to (MHD) that decay extremely fast. 

5.1 Proof of theorem 1.2 

Following [2], we define E as the space of all functions / G L^^^iM.'^) such that 

ll/IU = / \f{^)\dx + sup i? / \f{x)\dx (5.1) 

def7|a;|<l R>i J\x\>R 

is finite, and 

Jim R I \f\x)\dx = Q. 

Holder inequality implies that : 



^^+°° J\x\>R 



■ ^0 H >d+l (po < +oo) or 

Lll{W^)(ZE whenever <j Po 

'do>d + l (po = +oo). 

Let us prove that ||m||£; cannot remain uniformly bounded during a positive time interval, 
unless the orthogonality relations (1.9) are satisfied. 

Proposition 5.1 Let {u,B) e C{[0,T]; L'^{R'^) x ^^(R'^)) a solution to (MHD) such that 
uq E E. Assume that 

MeL°°([0,T];E) (5.2a) 

|M|2 + |5p gL°^([0,T];^). (5.2b) 
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Then there exists a constant c > such that the components of the initial data satisfy 

Vj, A; G {1, . . . , d}, [ ui4 - BiB', = c5j,k, (5.3) 

where Sj^k = 1 if j = k and otherwise. 

Proof. The proof will only be sketched briefly since it is a straightforward adaptation of [2]. 
Let us write the flrst equation of (MHD) in the following form (recall that S and Re can be 
set equal to 1): 



u{t) - e*^Mo 



+ V /" e(*-^)^ dj{u^u - B^B) ds = - f e(*-")^ VPis) ds, (5.4) 

•^1 Jo Jo 



where P = p + '-^ is the total pressure. Arguing as in [2], we see that (5.2) imply that all 
the terms in the left-hand side of (5.4) belong to L°°{[0,T]; E). Thus, we have 

VPgL°°([0,T];E) with P{t)= f e^'-'^^P{s)ds. 

Jo 

Let 

u^'^(t) = [ e^'^''^^u^u\s)ds 



Jo 
and 

&'\t)= I e^'-'^^B^B^{s)ds. 
Jo 

Taking the divergence in (5.4) yields 

d 

-AP = J2 9jdk{P''' - &'^). 
j,k=i 

One now deduces (5.3), applying Lemma 2.3 and Proposition 2.4 of [2]. D 

The proof of Theorem 1.2 is now very easy. Thanks to (1.3) and (1.4), assumptions (1.8a) 
and (1.8b) imply the existence of e' > e" > such that 

sup \u{t, ■)! e Ll C L}+,„ C E. 

te[o,T] 2+^+^ 

Moreover, the deflnition of the L^ decay rate at infinity (1.1) implies that 



and 



lim R''+^+^'' / \u{t,x)\^dx = 

■^^°° Jr<\x\<2R 



lim R^+'' / \B{t,x)\'^dx = 0, 

^^°° Jr<\x\<2R 



uniformly for t G [0,T]. Therefore 



sup (|M(t,-)|' + |i?(t,-)l') eL\^,„(lE. 

te[o,T] 



Conclusion (1.9) now follows from proposition 5.1. 
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5.2 Solutions of (MHD) with an exceptional spatial behavior 

We finally observe that solutions that decay faster than predicted by Theorem 1.3 do exist. 

Such solutions can be constructed starting with properly symmetric initial data. Assume, 
e.g., that Mq and Bq are rapidly decreasing in the usual pointwise sense when \x\ -^ +00 
(faster than any inverse polynomial) and that Auo{x) = Uo{Ax), ABq{x) = Bq{Ax) for all 
a; G M'^ and all matrix A E G, where G is a subgroup of the orthogonal group 0{d). Then 
the solution of (MHD) will inherit this property as far as it exists, the system being invariant 
under rotations. If the group G is rich enough, then these symmetry relations ensure the 
validity of conditions (1.9). Moreover the decay rate of the velocity field of the corresponding 
solution will depend on the symmetry group to which {uq, Bq) belongs. 

In dimension d = 2, 3 and for the Navier-Stokes equations, the optimal decay rates of the 
solution have been computed in [1] for each symmetry group. With simple modifications in 
the proofs, one could show that the same decay rates hold for the solution of (MHD). This 
is not surprising: indeed, since the magnetic field decays fast when \x\ -^ +00, the decay of 
the velocity field is governed only by the decay rate of the kernels Fj^^, defined by (2.1), and 
by the possible corresponding cancellations. These kernels are the same ones that appear in 
the Navier-Stokes system as well. 

Thus, for example, in dimension d = 2 and when G is the cyclic group of order n, one 
has 

Vt G [0,T*), u{t,x) = C(|a;|-("+i)) 

in the usual pointwise sense, when |x| -^ +00. In particular, the property of being simul- 
taneously completely invariant under rotations {i.e. G = SO {2)) and rapidly decreasing at 
infinity will be conserved by (m, B) during the evolution, if such property already holds for 
{uo,Bo). 

In dimension three, the largest decay rates of the velocity field {i.e. like |a;|~^ as x ^ +00) 
are obtained with the symmetry groups of the icosahedron. Those symmetric solutions are 
however unstable: in general, the velocity field of an infinitesimal perturbation of a highly 
symmetric fiow will decay much more slowly at infinity. 
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